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QGSJET-II: results for extensive air showers 
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The new hadronic Monte Carlo model QGSJET-II is applied for extensive air shower (EAS) calculations. The 
obtained results are compared to the predictions of the original QGSJET and of the SIBYLL 2.1 interaction 
models. It is shown that non-linear effects change substantially model predictions for hadron- nucleus interactions 
and produce observable effects for calculated EAS characteristics. Finally the impact of the new model on the 
interpretation of air shower array data is discussed. 



1. INTRODUCTION 

Experimental studies of high energy cosmic 
rays (CR) are mainly performed on the basis 
of extensive air shower (EAS) techniques: mea- 
suring various characteristics of nuclear-electro- 
magnetic cascades, induced by primary CR par- 
ticles in the atmosphere, one uses the obtained 
knowledge to reconstruct the properties of the 
initial particles. Because of a very complicated 
structure of the atmospheric cascade such an 
analysis necessarily requires both particle de- 
tectors - to measure various EAS components 
at Earth or in space, and proper simulation 
tools, which allow to mimic the cascade pro- 
cess and to establish a connection between the 
measured information and the primary particle 
characteristics. An important part of EAS sim- 
ulation programs are Monte Carlo (MC) mod- 
els of hadronic interactions. In the absence of 
a rigorous theoretical approach for a treatment 
of general (minimum-bias) hadron-nucleus and 
nucleus-nucleus collisions a development of reli- 
able hadronic MC generators for very high en- 
ergy interactions is far non-trivial; also due to 
the necessity to extrapolate corresponding exper- 
imental knowledge from the energies of present 
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colliders to much higher CR energies. On the 
other hand, measuring different EAS components 
allows to perform a cross-check of model descrip- 
tion and to discriminate between available MC 
generators. 

During last years QGSJET model [T05] has 
been successfully used in the field, in particu- 
lar, being employed in popular EAS simulation 
programs, CORSIKA |1| and AIRES 0. Also 
model tests, performed for example by the KAS- 
CADE collaboration, using its multi-component 
EAS detector setup [H], showed that the model 
reproduces measured shower characteristics and 
their correlations reasonably well |7l8j . However, 
the mentioned agreement is still far from being 
ideal; recent investigations showed some discrep- 
ancies in the description of basic correlations be- 
tween electron and muon EAS components [5]. 
Also annoying is the fact that analyzing different 
EAS components results in sizably divergent con- 
clusions concerning primary CR composition, for 
all models in use |1UI11| . 

Recently a new hadronic MC model QGSJET- 
II has been developed ^21; which included a treat- 
ment of non-linear interaction effects in individ- 
ual hadronic (nuclear) collisions. The latter gave 
a possibility to obtain a consistent description 
of interaction cross sections and parton momen- 
tum distributions in hadrons, compared to cor- 
responding measurements, thus solving the basic 
deficiency of the original QGSJET and providing 
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a much more reliable model extrapolation into 
ultra-high energy domain. Additionally, realis- 
tic nuclear density parameterizations have been 
employed in the model, individually for each nu- 
cleus |13| . more reliable two-component low mass 
diffraction treatment has been used and all 
the model parameters have been re-calibrated us- 
ing a wider set of accelerator data. The goal of 
the present work is to investigate the impact of 
the mentioned modifications on the calculated air 
shower characteristics and to draw possible con- 
sequences for EAS data interpretation. 




Figure 1. Inelastic hadron-air cross sections for 
QGSJET-II, QGSJET-I, and SIBYLL 2.1 models 
- smooth, dashed, and dotted curves. 



2. RESULTS FOR HADRON-NUCLEUS 
INTERACTIONS 

Despite the fact that the physics of hadronic 
interactions is very complicated and contributes 
to air shower development in a rather non-trivial 
way, basic EAS observables mainly depend on 
a limited number of macroscopic characteristics 
of hadron-air collisions. Indeed, the most fun- 
damental EAS parameter, the position of the 
shower maximum ATmax is grossly defined by in- 
elastic hadron-air cross sections and by the in- 
teraction inelasticities, the latter quantity being 
defined as the relative difference between the lab. 
energies of the initial and the most energetic fi- 
nal particles. While the number of charged lep- 
tons Ne, measured at a given observation level, is 
strongly correlated with X^axj this correlation is 
much weaker for corresponding number of muons 
N^, and even becomes negative for relatively high 
muon threshold. On the other hand, muon num- 
bers appear to be quite sensitive to multiplicities 
of elementary interactions. 

In Figs, no the predictions of the new model 
for inelastic cross sections a™^\^j., inelasticities 
-ft'jjifairi ^^'^ multiplicities of charged particles 
-^h^ain ^^''^ plotted in comparison with both the 
original QGSJET and the SIBYLL 2.1 [HI mod- 
els. Compared to the original model, the re- 
sults of the QGSJET-II are defined by a com- 
petition of two effects: non-linear screening cor- 
rections, which lead to a reduction of the interac- 
tion eikonal and correspondingly of the number 
of elementary particle production processes (cut 
Pomerons), and steeper parton momentum dis- 




Figure 2. Inelasticities of hadron-air interac- 
tions for QGSJET-II, QGSJET-I, and SIBYLL 
2.1 models - smooth, dashed, and dotted curves. 



tributions, leading to a faster energy increase of 
the latter quantities ^^l- At not too high ener- 
gies the first effect dominates, leading to smaller 
numbers of cr™,!:^!!. , K}^^]^-^., and A^^'iair> especially 
for pion-air interactions. On the other hand, in 
the very high energy limit the infiuence of parton 
distributions prevails and the new model predicts 
larger values for the quantities of interest. 

3. AIR SHOWER CHARACTERISTICS 

Sizably smaller inelasticities of the new model 
lead to a somewhat deeper position of the shower 
maximum, compared to the original QGSJET - 
Fig. ^ The corresponding shift of X^^ax changes 
from about 20 g/cm^ at lO^'^ eV to ~10 g/cm^ at 
10^*^ eV in case of proton-induced EAS; for iron- 
induced EAS it increases slowly from 15 g/cm^ till 
~20 g/cm^. As the relative strength of non- linear 
effects is larger for nucleus-nucleus collisions one 
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Figure 3. Multiplicities of charged parti- 
cles in hadron-air interactions for QGSJET-II, 
QGSJET-I, and SIBYLL 2.1 models - smooth, 
dashed, and dotted curves. 
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Figure 5. Relative iVe-difference at sea level be- 
tween QGSJET-II and QGSJET-I models for ver- 
tical proton- and iron-induced EAS - smooth and 
dashed curves. 
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Figure 4. Average position of the shower max- 
imum for proton- and iron-induced EAS for 
QGSJET-II, QGSJET-I, and SIBYLL 2.1 models 
- smooth, dashed, and dotted curves. The points 
are the data of the HIRES collaboration |15j . 
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Figure 6. Relative A'^f^-difference {E^, > 10 GeV) 
at sea level between QGSJET-II and QGSJET- 
I models for vertical proton- and iron-induced 
EAS - smooth and dashed curves, and be- 
tween SIBYLL 2.1 and QGSJET-I models for 
p— induced EAS - dotted curve. 



observes sizable deviations from the superposition 
picture for nucleus-induced air showers. 

The relative difference between QGSJET-II 
and QGSJET-I models for predicted numbers of 
electrons and muons (i?^ > 10 GeV) at sea level 
for vertical proton- and iron-induced showers is 
shown in Figs. El El for comparison the differ- 
ence between SIBYLL 2.1 and QGSJET-I is also 
shown for the case of primary protons. While for 
electron numbers we observe the expected corre- 
lation with the position of the shower maximum, 
the predicted muon numbers in the new model 
are significantly smaller, by as much as 30% at 
highest energies, due to a strong reduction of in- 
teraction multiplicity, especially for pion-nucleus 



interactions, with the results being close to the 
predictions of the SIBYLL 2.1 model. 

4. DISCUSSION 

Systematic account for non-linear screening ef- 
fects in the QGSJET-II model results in a sub- 
stantial reduction of multiplicity and inelasticity 
of hadron-nucleus interactions compared to the 
original QGSJET. This leads in turn to a shift 
of the predicted position of the shower maxi- 
mum deeper in the atmosphere, to larger num- 
ber of electrons and smaller number of muons 
at sea level. Applying the new model to EAS 
data reconstruction will change present conclu- 
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sions concerning CR composition towards heav- 
ier primaries. The corresponding X^ax change 
is not large at highest energies. Nevertheless, 
while the predictions of the original QGSJET, 
being compared to HIRES measurements, are 
marginally consistent with the assumption of ul- 
tra high energy cosmic rays being only protons, 
this no longer the case with the new model. 

Similarly, at the knee energies QGSJET-II pre- 
dicts a substantially smaller number of muons for 
a given number of electrons, which would lead 
to a conclusion about a sizably heavier composi- 
tion, compared to what is currently obtained with 
QGSJET-I. In particular, the mentioned change 
should bring closer together the composition re- 
sults of the KASCADE collaboration obtained 
on the basis of electron/muon and hadronic EAS 
components ^01 ■ Also the obtained changes seem 
to go in the right direction to resolve the reported 
discrepancies in the electron/muon correlations 
[S]. It should be mentioned, however, that at 
highest energies the disagreement of the composi- 
tion results obtained with fluorescence light-based 
measurements and with ground arrays ^l] would 
increase using the new model, as the reduction 
of EAS muon number at highest energies is much 
stronger compared to the corresponding effect for 
the shower maximum. 

It is of interest to investigate the impact of the 
new model on the experimental techniques of pri- 
mary energy reconstruction, where two methods 
are mainly used: either based on the integral of 
the shower cascade curve in case of fluorescence 
light-based measurements or estimating the en- 
ergy from charged particle densities at distances 
600 and 1000 m from the shower core JH| ■ In the 
first case model dependence of the results enters 
mainly via predicted numbers of charged particles 
in the shower maximum N^^'^, which shows very 
small model sensitivity - less than 3% between 
QGSJET-II, QGSJET-I, and SIBYLL 2.1. Simi- 
larly, for electron density at large distances from 
the core one observes just a few per cent shift be- 
tween the new and the old versions of QGSJET. 
However, depending on experimental techniques 
applied, a substantial contribution to the mea- 
sured signal may come from the muon compo- 
nent, which is the case for example, for AUGER 



collaboration T7|. As the latter is substantially 
reduced in QGSJET-II, corresponding energy es- 
timates of the ground array would move upwards 
compared to QGSJET-I. On the other hand, the 
possibility to perform independent energy recon- 
structions with both the fluorescence detector and 
the ground array could allow to perform a model 
consistency check at the highest CR energies. 
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